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Abstract — In heterogeneous cellular networks (HCNs), it is 
desirable to offload mobile users to small cells, which are typically 
significantly less congested than the macrocells. To achieve suffi- 
cient load balancing, the offloaded users often have much lower 
SINR than they would on the macrocell. This SINR degradation 
can be partially alleviated through interference avoidance, for 
example time or frequency resource partitioning, whereby the 
macrocell turns off in some fraction of such resources. Naturally, 
the optimal offloading strategy is tightly coupled with resource 
partitioning; the optimal amount of which in turn depends on 
how many users have been offloaded. In this paper, we propose 
a general and tractable framework for modeling and analyzing 
joint resource partitioning and offloading in a two-tier cellular 
network. With it, we are able to derive the downlink rate 
distribution over the entire network, and an optimal strategy 
for joint resource partitioning and offloading. We show that load 
balancing, by itself, is insufficient, and resource partitioning is 
required in conjunction with offloading to improve the rate of 
cell edge users in co-channel heterogeneous networks. 



I. Introduction 

The exponential growth in mobile traffic primarily driven by 
mobile video [1] has led the drive to deploy low power base 
stations (BSs) both in licensed and unlicensed spectrum in 
order to complement the existing macro cellular architecture. 
Increased cell density increases the area spectral efficiency 
of the network |2| and is the key factor for the very large re- 
quired capacity boost. Such a heterogeneous network (HetNet) 
consists of macro BSs coexisting with small cells formed by 
co-channel low power base stations like micro, pico and femto 
BSs, as well as unlicensed band WiFi access points (APs) |3), 

The "natural" association/coverage areas of the low power 
APs tend to be much smaller than those of the macro BSs, 
and hence the fraction of user population that is offloaded to 
small cells may often be limited, resulting in insufficient relief 
to the congested macro tier. This user load disparity not only 
leads to suboptimal rate distribution across the network, but the 
lightly loaded small cells may also lead to bursty interference 
causing QoS degradation (5), One practical technique 
for proactively offloading users to small cells is called cell 
range expansion (CRE) [3| wherein the users are offloaded 
through an association bias. A positive association bias implies 
that a user would be offloaded to a small cell as soon as 
the received power difference from the macro and small cell 
drops below the bias value, which "artificially" expands the 
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association areas of small cells^J Though experiencing reduced 
congestion, in co-channel deployments such offloaded users 
also have degraded signal-to-interference-ratio (SINR), as the 
strongest AP (in terms of received power) now contributes to 
interference. Therefore, the gains from balancing load could 
be negated if suitable interference avoidance strategies are not 
adopted in conjunction with cell range expansion particularly 
in co-channel deployments [6]. One such strategy of interfer- 
ence avoidance is resource partitioning |3), (7), wherein the 
transmission of macro tier is periodically muted on certain 
fraction of radio resources (also called almost blank subframes 
in 3GPP LTE). The offloaded users can then be scheduled in 
these resources by the small cells leading to their protection 
from co-channel macro tier interference. 

A. Motivation and Related Work 

It has been established that without proactive offloading 
and resource partitioning only limited performance gains can 
be achieved from the deployment of small cells j8)-fl2). 
These techniques are strongly coupled and directly influence 
the rate of users, but the fundamentals of jointly optimizing 
offloading and resource partitioning are not well understood. 
For example, an excessively large association bias can cause 
the small cells to be overly congested with users of poor 
SINR, which requires excessive muting by the macro cell 
to improve the rate of offloaded users. Earlier simulation 
based studies 1 1 1 1, p2) confirmed this insight and showed 
that excessive biasing and resource partitioning can actually 
degrade the overall rate distribution, whereas the choice of 
optimal parameters can yield about 2-3x gain in the rate 
coverage (fraction of user population receiving rate greater 
than a threshold). Although encouraging, a general tractable 
framework for characterizing the optimal operating regions for 
resource partitioning and offloading is still an open problem. 
The work in this paper is aimed to bridge this gap. 

A "straightforward" approach of finding the optimal strat- 
egy is to search over all possible user-AP associations and 
time/frequency allocations for each network configuration. 
Besides being computationally daunting, this approach is 
unlikely to lead to insight into the role of key parameters 
on system performance. Another methodology is a probabilis- 
tic analytical approach, where the network configuration is 
assumed random and following a certain distribution. This 
has the advantage of leading to insights on the impact of 
various system parameters on the average performance through 

'Access point (AP), base station (BS), cell are used interchangeably in the 
paper. 
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tractable expressions. Analytical approaches to study biasing 
and interference coordination were studied in (T3j-|[l5j, but 
downlink rate (one of the key metrics) was not investigated. 
Optimal bias and almost blank subframes were prescribed in 
p3) based on average per user spectral efficiency. A related 
SINR based analysis was done in fl4) for resource partitioning 
in frequency domain, but offloading was not captured. The 
choice of optimal range expansion biases in fT3| was not based 
on rate distribution. In this paper, we use the metric of rate 
coverage, which captures the effect of both SINR and load 
distribution across the network. Semi-analytical approaches in 
[16 1, 1 17 1 showed, through simulations, that there exists an 
optimal association bias for fifth percentile and median rate 
which is confirmed in this paper through our analysis. Also, to 
the best of our knowledge, none of the mentioned earlier works 
considered the impact of backhaul capacities on offloading, 
which is another contribution of the presented work. 

B. Approach and Contributions 

We propose a general and tractable framework to analyze 
joint resource partitioning and offloading in a two-tier cellular 
network in Section]!!] The proposed modeling can be extended 
to a multiple tier setting. Each tier of base stations is mod- 
eled as an independent Poisson point process (PPP), where 
each tier differs in transmit power, path loss exponent, and 
deployment density. The mobile user locations are modeled 
as an independent PPP and user association is assumed to 
be based on biased received power. On all channels, i.i.d. 
Rayleigh fading is assumed. Similar tractable frameworks were 
used for deriving SINR distribution in HCNs in [18|-|20). The 
empirical validation in [21] and theoretical validation in [22| 
for heavily shadowed cellular networks have strengthened the 
case of modeling macro cellular networks using a PPP. Due 
to the formation of random association/coverage areas in such 
network models, load distribution is difficult to characterize. 
An approximate load and rate distribution was derived for 
multiple radio access technology (RAT) HetNets in (23]. 

Based on our proposed approach, the contributions of the 
paper can be divided into two categories: 
Analysis. The rate complementary cumulative distribution 
function (CCDF) in a two-tier co-channel heterogeneous 
network is derived as a function of the cell range expan- 
sion/offloading and resource partitioning parameters in Section 



TABLE I: Notation Summary 



III Rate coverage at a particular rate threshold is the rate 
CCDF value at that threshold. The derived rate distribution 
is then modified to incorporate a network setting where APs 
are equipped with limited capacity backhaul. Under certain 
plausible scenarios, the derived expressions are in closed form. 
Design Guidelines. The theoretical results lead to joint re- 
source partitioning and offloading insights for optimal SINR 
and rate coverage in Section [IV] In particular, we show the 
following: 

• With no resource partitioning, optimal association bias 
for rate coverage is independent of the density of the 
small cells. In contrast, offloading is shown to be strictly 
suboptimal for SINR in this case. 

• With resource partitioning, optimal association bias de- 
creases with increasing density of the small cells. 



Notation 


Description 


<I> fe ; § u 


PPP of APs of fc th tier; PPP of mobile users 




Density of APs of k th tier; density of mobile users 


r k) ^k 


Transmit power of APs of k^ tier; normalized transmit 
power of APs of fc th tier 


Bfc; B fe 


A • !_■ J? 1 ||l i' 1* J ' 4-' 1_" 

Association bias tor k tier; normalized association bias 
for k th tier. 


ctk\ &k 


Path loss exponent of fc th tier; normalized path loss 
exponent of k' h tier 


W;O fe 


Air interface bandwidth at an AP for resource allocation; 
backhaul bandwidth at an AP of fc th tier 


Ui 


Macro cell users 1 = 1, small cell users (non-range 
expanded) / = B, offloaded users I = B 


mn 


Resource partitioning fraction; inverse of the effective 
fraction of resources available for users in lAi 




Map from user set index to serving tier index, J(l) = 1, 


■2 

u 


Thermal noise power 


■Ai 


Association probability of a typical user to XAi 


1l;S;p 


Rate coverage; SINR coverage; rate threshold 




Load at tagged AP of u 6 W; ; PMF of load 

p l (n)=P(N l = n) 


Z k ;Y t 


Distance of the nearest AP in k th tier; distance of the 
tagged AP conditioned on u £ Ui 


Cx k ; C'k 


Association region; area of an AP of tier k 



• In both of the above scenarios, the optimal fraction 
of users offloaded, however, increases with increasing 
density of small cells. 

• With decrease in backhaul capacity/bandwidth the opti- 
mal association bias for the corresponding tier always 
decreases. However, in contrast to the trend in the "in- 
finite']^] backhaul scenario, the optimal association bias 
may increase with increasing small cell density. 

The paper is concluded in Section [V] and future work is 
suggested. 

II. Downlink System Model and Key Metrics 

In this paper, the wireless network consists of a two-tier 
deployment of APs. The location of the APs of k th tier (k = 
1,2) is modeled as a two-dimensional homogeneous PPP 
of density (intensity) Afc. Without any loss of generality, let the 
macro tier be tier 1 and the small cells constitute tier 2. The 
locations of users (denoted by U) in the network are modeled 
as another independent homogeneous PPP $ u with density \ u . 
Every AP of fc th tier transmits with the same transmit power 
Pfe over bandwidth W. The downlink desired and interference 
signals from an AP of tier-fc are assumed to experience path 
loss with a path loss exponent a^. A user receives a power 
PkH x x~° k from an AP of k th tier at a distance x, where 
H x is the random channel power gain. The random channel 
gains are assumed to be Rayleigh distributed with average 
unit power, i.e., H x ~ exp(l). General fading distributions 
can be considered at some loss of tractability ]24) . The noise 
is assumed additive with power a 2 . The notations used in this 
paper are summarized in Table [I] 

2 Infinite bandwidth implies sufficiently large so as not to affect the effective 
end-to-end rate. 
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A. User Association 

The analysis in this paper is done for a typical user u located 
at the origin. This is allowed by Slivnyak's theorem |25|, 
which states that the properties observed by a typicaj^] point of 
a PPP $ is same as those observed by a node at origin in the 
process $U{0}. Let Z k denote the distance of the typical user 
from the nearest AP of fc th tier. It is assumed that each user 
uses biased received power association in which it associates 
to the nearest AP of tier j if 



j = arg max P k B k Z, ak , 
fce{i,2} 



(1) 



where is the association bias for k tb tier. Increasing 
association bias leads to the range expansion for the corre- 
sponding APs and therefore offloading of more users to the 
corresponding tier. For clarity, we define the normalized value 
of a parameter of a tier as its value divided by the value it 
takes for the serving tier. Thus, 

a a Pfe a a B fc „ A a k 

P/c = — , B fc = — , and a k = — , 
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are respectively the normalized transmit power, association 
bias, and path loss exponent of tier k conditioned on the user 
being associated with tier j. Association bias for tier 1 can be 
assumed to be unity (Bi = dB) and that of tier 2 be simply 
denoted by B. In the given setup, a user u G U can lie in 
following three disjoint sets: 

{Ui if j = 1, PiZ~ ai > P 2 B^ a! 
U B if j = 2 and P 2 Z^ a2 > PiZf ai 
U B if ? = 2 and P 2 Z^ a2 < PiZf ai < P 2 BZ^ a2 , 

(2) 

where U\ U Ub U U§ = U clearly. The set U\ is the set of 
macro cell users and the set Ub is the set of unbiased small 
cell users. Thus, the set Ub is independent of the association 
bias. The users offloaded from macro cells to small cells due 
to cell range expansion constitute Ub and are referred to as the 
range expanded users. All the users associated with small cells 
are U2 = U§ U Ub- We define a mapping J : {1, B, B} —> 
(1, 2} from user set index to serving tier index. Thus, from 
(gj, J(l) = 1, J(B) = J(B) = 2. 

The biased received power based association model de- 
scribed above leads to the formation of association/coverage 
areas in the Euclidean plane as described below. 

Definition 1. Association Region: The region of the Eu- 
clidean plane in which all users are served by an AP is called 
its association region. Mathematically, the association region 
of an AP of tier j located at x is 



C Xj = \y G R 2 : \\y-x\\ < 



where XZ(y) — arg min \\y 



PfeB fe 



\y-X* k (y)\\ & *V kj 
(3) 



The random tessellation formed by the collection {C Xj } of 
association regions is a general case of the circular Dirich- 
let tessellation [26 1. The circular Dirichlet tessellation (also 



known as multiplicatively weighted Voronoi |27| Chapter 3]) 
is the special case of the presented model with equal path loss 
coefficients. 

B. Resource Partitioning 

A resource partitioning approach is considered in which 
the macro cell shuts its transmission on certain fraction of 
time/frequency resources and the small cell schedules the 
range expanded users on the corresponding resources, which 
protects them from macro cell interference. 

Definition 2. rj: The resource partitioning fraction 77 is the 
fraction of resources on which the macro cell is active, where 

< T] < 1. 

Thus, with resource partitioning 77 fraction of the resources 
at macro cell are allocated to users in U\ and those at small 
cell are allocated to users in Uq. The fraction 1 — 77 of the 
resources in which the macro cell shuts down the transmission, 
the small cells schedules the range expanded users, i.e., Ub- 
Let 7/ denote the inverse of the effective fraction of resources 
available for users in Ui. Then, 7; = 1/?/ for I £ {l,B} 
and 7j = 1/(1 — 77) for I = B. The operation of range 
expansion and resource partitioning in a two-tier setup is 
further elucidated in Fig. [T] In these plots, the power ratio 

is assumed to be &l = 20 dB and B = 10 dB. 
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As a result of resource partitioning (0 < 77 < 1), the SINR 
of a typical user u, when it belongs to Ui, is 



sim = i(le{i,B } ) ^ H f\ 
E fe =i h + v 2 



-1(1 = B) 



P 2 H y y^ 
/2+a 2 ' 
(4) 

where 1(A) denotes the indicator of the event A, H y is the 
channel power gain from the tagged AP s/ (AP serving the 
typical user) at a distance y, Ik denotes the interference from 
the fc th tier. The set of APs contributing to interference from 
the k th tier are from $/. \ s^. Thus 



h = p k H * x ~ ak 

x£<S> k \si 



(5) 



Proportional rate allocation is assumed at each AP where 
the rate allocated to each associated user is proportional 
to its spectral efficiency. Further, user queues are assumed 
saturated implying that each AP always has data to transmit 
to its associated mobile users. Let U s denote the set of users 
associated with the tagged AP. If the tagged AP belongs to 
macro tier, then N\ — \U S n U\ \ denotes the total number of 
users (or load henceforth) sharing the available r\ fraction of 
the resources. Otherwise, if the tagged AP belongs to tier 2, 
then the load is N 2 = \U s nU 2 \ of which N B = \U s f~)Ug\ users 
share the r\ fraction of the resources and Nb = \U s DUb \ users 
share the rest 1 — 77; N 2 = Nb + N§ — 1 (one is subtracted 
to account for double counting of the typical user). Thus, the 
rate of a typical user u is 



R = { l(u G Ui 



The term typical and random are interchangeably used in this paper. 



i(ueU B )^+l(ueU B 
x Wlog(l 



1 



n 



N B 
SINR) . 



(6) 



(a) Active macro tier 



(b) Muted macro tier 



Fig. 1: A filled marker is used for a node engaged in active transmission (BS) or reception (user), (a) The macro cells (filled squares) serve 
the macro users U\ and small cells (filled diamonds) serve the non-range expanded users (Wg) (filled circles), (b) The macro cells (hollow 
square) are muted while the small cells (filled diamonds) serve the range expanded users Ub (filled circles in the shaded region). 



The above rate allocation model assumes infinite backhaul 
bandwidth for all APs, which may be particularly questionable 
for small cells. Discussion about limited backhaul bandwidth 
is deferred to Sec. IIII-DI 

C. Rate and SINR Coverage 

The rate and SINR coverage can be formally defined as 
follows. 

Definition 3. Rate/SINR Coverage: The rate coverage for a 
rate threshold p is 

K(p)±V(R>p), (7) 



and SINR coverage for a threshold t is 
S(r) = P(SINR > t) 



(8) 



The coverage can equivalently interpreted as (i) the prob- 
ability that a randomly chosen user can achieve a target 
threshold, (ii) the average fraction of users in the network who 
at any time achieve the corresponding threshold, or (iii) the 
average fraction of the network area that is receiving rate/SINR 
greater than the rate/SINR threshold. 

III. Rate Distribution 

This section derives the rate distribution (coverage) and is 
the main technical section of the paper. 

A. Association Areas and Load Distribution 

For completely characterizing the rate distribution, the aver- 
age fraction of users belonging to the respective three disjoint 
sets QA\, Ujj, and Ub) is needed. Using the ergodicity of the 
PPP, these fractions are equal to the association probability of 
a typical user to these sets, which are derived in the following 
lemma. 



Lemma 1. (Association probabilities) The association prob- 
ability, defined as Ai = P(it G U\), is given below for each 
set 



)2/a h ^/& k d 



A 1 =2ir\ 1 zcxp -7T V Afc(P fc B fe ) ; 
Jo \ k=i 

A B =2Tr\ 2 J^ zexp (-tt £ ^k(Pk) 2/ak z 2/&k ^j dz. 



(9) 



(10) 

A B = 27r\ 2 j^ zjexp (-n ^ A fe (P fc B fe ) 2/a ^ 2/< ^ 

~ ex P [~KY,M?k) 2/ak z 2/& ^j jdz. (11) 

If path loss coefficients are same, i.e., = a, the association 
probabilities simplify to: 



A, 



Ar 



Ai 



,A T 



A 2 



ELiA fe (P fc B fc ) 2 /« ' B £LiA fc (P fc ) 2 A* 



E.=iA fe (P fc B fc )V« ELiA fc (P fe ) 2 / Q 



(12) 



Proof: See Appendix |A] ■ 
Equation ( fT2| ) corroborates the intuition that increasing 
association bias B leads to decrease in the mean population of 
macro cell users implied by the decreasing A\. On the other 
hand, the mean population of range expanded users increases 
implied by the increasing A 2 °- Further, Ai — A§ + Ab is 
the probability of a typical user associating with the tier 2. 

Remark 1. {Mean Association Area) Association area of an 
AP is the area of the corresponding association region. Using 
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the ergodicity of the PPP, the mean of the association area 
C k of a typical AP of k th tier is E [C k ] = 

The association region of a tier 2 AP can be further 
partitioned into two regions. The non-shaded region in Fig. 
[T] surrounding a small cell at x can be characterized as 

C XB ± {y e K 2 : \\y-x\\ < (P 2 /P 1 ) 1/a2 \\y-X* k (y)\\ & \ Vfc}. 

(13) 

As per (|2jl, all the users lying in C xs are the small cell users 
(belonging to Uq) and recalling ^ all users lying in C XB = 
C X2 ~ Cx s are me offloaded users that belong to Ub- In Fig. 
[T| C Xb is the shaded region surrounding a tier 2 AP. 

Remark 2. (Association Area Distribution) A linear scaling 
based approximation for the distribution of association areas 
proposed in l{23^ , which matched the first moment, is general- 
ized in this paper to the setting of resource partition as below 



d = C 



Ai 
Ax)' 



C U C { ^),andC B = c(^- 



(14) 
(15) 



where C (y) is the area of a typical cell of a Poisson Voronoi 
(PV) of density y (a scale parameter). 

Using the area distribution proposed in J28| for PV C(y), 
the following lemma characterizes the probability mass func- 
tion (PMF) of the load seen by a typical user. 

Lemma 2. (Load PMF) The PMF of the load at tagged AP 
of a typical user u *E Ui is 



p l (n)^¥(N l =n) 

3.5 3 5 r(n + 3.5) f\ u A 



(n-1)! T(3.5) V A 



V(i) 



3.5 



n > 1, (16) 



where T(x) = L exp(— t)t x 1 dt is the gamma function. 
Proof: See Appendix [B] ■ 

B. SINR Distribution 

The conditional SINR coverage, when a typical user u 
is Si{t) = P(SINR > t\u eUi) . 

Lemma 3. fSINR Coverage) For a typical user in the setup 
of Sec. [7/] the SINR coverage is 

S{t) = AiSi(t) + A b S b {t) + A b S b {t), (17) 

where the conditional SINR coverage are 



OO f 

S\(t) = 2n— - / uexp< 

W Ai J | SNRi 

n v 



(y) 



7r^A fe P^ /afc Q(r,a fel B fc y/^ \dy (18) 



SNR 2 (y) 



- 7T A fe Pfc /afc Q(r, a k , l)y 2 / & * \dy (19) 

k=l ) 



Sb(t 

— 7rAl(PlBl) 

Q(a, 6, c) 



2 ^y ^ exp 



2/ai 2/a 

2/ 



{ 

1 1 - exp | 



SNR 2 (y) 



7rA 2 Q(r, a 2 ,l)y 2 



„2/b 



smjy)-*^' *' 1)2/2 

TrAiCPi) 2 ^ 1 ^ 1 }!^, (20) 
and SNR fe (y) = 



(f ) 2/6 l+u b / 2 



Proof: See Appendix [C] ■ 
The result in Lemma [3] is for the most general case and 
involves a single numerical integration along with a lookup 
table for Q. The expressions can be further simplified as in 
the following corollary. 

Corollary 1. With noise ignored, SNR^, — > oo, and equal path 
loss coefficients a k = a, the SINR coverage of a typical user 



S(t) 



Ai 



ELi A fc (P fc /P!) 2 /"Q(r, a, B fc /Bi) 

+ h 

ELi^(P fc /P2) 2 ^Q(r,a,l) 

+ *3 

A 2 Q(r, a, 1) + Ai {PiB 1 /(P 2 B 2 )} 2/c 

A 2 



A 2 Q(r,a,l) + Ai(Pi/P 2 ) 2 /° 



(21) 



As evident from the above Lemma and Corollary, SINR 
coverage is independent of the resource partitioning fraction 
T) because of the independence of SINR on the amount of 
resources allocated to a user in our model. Further, the SINR 
distribution of the small cell users, S§, is independent of 
association bias, as Ub is independent of bias. Further insights 
about SINR coverage are deferred until the next section. 
In general, we show that SINR coverage with and without 
resource partitioning show considerably different behavior, 
which is also reflected in the rate coverage trends. 

C. Main Result 

Similar to the conditional SINR coverage, conditional 
rate coverage, when a typical user u £ Ui is IZi(p) = 
P (R > p\u e Ui) . The following theorem gives the rate dis- 
tribution over the entire network. 

Theorem 1. (Rate Coverage) For a typical user in the setup 
of Sec. [7/] the rate coverage is 

K{p) = AyRi ( P ) +A B K B {p) + AbTIb (p) , (22) 



k=l 



where the conditional rate coverage are given by ( |23| )-(25 1, 

t{x) = 2 X - 1, and p = p/W. 
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= 2tt^- J™ V exp | - - A fe P^ fe Q(«(np 7 i), «*, B fe )y 2/Sfc jdj, (23) 

= 27r-^^^(n) r y exp ( - -^^A fe P^ fc Q(^ 7i5 ),« fc ,l)y 2/,i 4d y (24) 

B n>l ^° I 2 ^ k=l ) 

K B (p) = 2^ 5> B (n) J exp {"|^gy - 7rA 2 y 2 Q(t(np 7s ), a 2 , 1) - ^A l2 ; 2 /^ (P 1 B 1 ) 2 ^ } 

- exp {-^g|y ~ 7r\ 2 y 2 Q(t(np lB ), a 2 , 1) - ■n\\y 2 ^ &1 (Px) 2 ^ } jdy, (25) 

^(p^TT^jTyexpj-^^ (29) 

fc fi (p) = 2^ ^ y exp | - ^fg 1 - * E A fc P' /Qfc Q(i(^57B), a fc , I)?/ 2 /- jdy (30) 
7l B (p) = 2^ J™ y j exp j-^ggl - n\ 2 y 2 Q(t(pN B ~/ B ), a 2 , 1) - ttAu/ 2 /" 1 (P 1 B 1 ) 2 /^ } 

~ 6XP Sy) } ~ * X *V 2 W(pXb'Tb), <*2, 1) - 7rA iy 2 /^ (P0 2/ai | U, (31) 



¥(R > p)= P(u G Wi)P I ^ log(l + SINR) > p|u € 
( /« •:- ///. i? ( v - log(l + sinr) > P \u e U B 

(u G W B )P I " v ~ '" log(l + SINR) > p\u G W B 



Proof: Using Q and (|6j, the probability that the rate 
requirement of a random user u is met is 

'W(l-Tj) 

(26) 

= [5i (t(pJViTi))] +^ fl E^ [S B (t(pN BlB ))} 

+ A B E NB [S B (t(pN BlB ))}, (27) 

where p = p/W. Using Lemma ^ E N , [Si(t(xNi))} = 
En>iB( n )^i (t(x?i)). Using the law of total probability and 
Lemma [3] the rate coverage expression is obtained. ■ 

The rate distribution expression for the most general set- 
ting requires a single numerical integral after use of lookup 
tables for Q and V. The summation over n in Theorem [T] 
can be accurately approximated as a finite summation to a 
sufficiently large value, n max (say), since both the terms pi (n) 
and Si (t(xn))) decay rapidly for large n. 

The rate coverage expression can be further simplified if the 
load at each AP is assumed to equal its mean. 

Corollary 2. (Mean Load Approximation) Rate coverage with 
the mean load approximation is given by 

TZ(p) = AiHi{p) + A B TZ B {p) + A B K B (p), (28) 

where the conditional rate coverage are given by \29\ -\31\ 

and Ni=E \NA = 1 + 1 - 28A "- 4 ' . 



Proof: Lemma |2] gives the first moment of load as 
E [Ni] = 1 + ^ME[C 2 (1)]. Further, using the result 

that E [C 2 (l)j = 1.28 0291, along with an approximation 
E Nk [S k (t(xN k ))] « Sklt(xE[N k })), the simplified rate 
coverage expression is obtained. ■ 
The mean load approximation above simplifies the rate 
coverage expression by eliminating the summation over n. The 
numerical integral can also be eliminated by ignoring noise 
and assuming equal path loss exponents (as is done in Sec 
|IV-B[ ). As can be observed from Theorem [T] and Corollary [2] 
the rate coverage for range expanded users 1Z B increases with 
decrease in resource partitioning fraction rj, as users in U B 
can be scheduled on a larger fraction of (macro) interference 
free resources. On the other hand, the rate coverage for the 
macro users 1Z 1 and small cell (non-range expanded) users 7Z B 
decreases with the corresponding decrease. Further insights on 
the effect of biasing are delegated to the next section. 

D. Rate Coverage with Limited Backhaul Capacities 

Analysis in the previous sections assumed infinite backhaul 
capacities and thus the air interface was the only bottleneck 
affecting downlink rate. However, with limited backhaul ca- 
pacities 0^ for BSs of tier k, the rate is given by 

R' = l(ue Wi)min ( R, ) + l(u G U 2 )min ( R, 

(32) 

where R is the rate of the user with infinite backhaul band- 
width. The above rate allocation assumes that the available 
backhaul bandwidth for a BS of tier k, Ofe, is shared equally 



7 



among the associated users/load Nk- The following lemma 
gives the rate distribution in this setting. 

Lemma 4. (Rate Coverage with Limited Backhaul) The rate 
coverage in the setting of Sec. ^and with rate model of {32 > 



n'(p) =F(R' >p) = A^ip) + A B 1l B {p) + A B 1l B {p), 

(33) 

where 

roi/p-ii 

n 'x(p)= Y Pi( n ) s i (*(7i«P)) ) 

n=l 

[0 2 /p~2] [0 2 /p-rn-l] 

n s(P)= Y PB{m) P§{ n )S§ (ths n P)) > 

m— n— 1 

ro 2 /p-2] rOa/p-m-ll 

1Z 'b(p)= Y Ps( m ) Y PB{n)S B {t( lB np)) , 

m=0 n=l 

(34) 

and Si is given by Lemma [3] 

Proof: Since the maximum rate of a user u € lii is 
Ojn\/Njm. Thus, for this user to have positive rate coverage, 
i.e., P(R > p) > 0, a necessary condition is iV/m < 



r- 



1]. When this necessary condition is satisfied, the 



rate coverage is equivalent to 

Til' = P (r > p) 



^ log(l + SINR) > p n N J{1) < ^ ] (35 I 



li N, 



P 



Using Ni — Njn) for I = 1, and independence of N B and 
N B the conditional rate coverage are obtained. ■ 

It is evident from the above Lemma that rate coverage 
decreases with decreasing backhaul bandwidth. Therefore, 
decreasing O2 will lead to decrease in the rate of the user 
when it is associated to small cell and thus decreasing the 
optimal offloading bias (this is further explored in subsequent 
sections). As the backhaul bandwidth increases to infinity, 
Lemma Q leads to Theorem [lj or, limo J(!) ->-oo f^i — > T^i- 

E. Validation of Analysis 

We verify the developed analysis, in particular Theorem [TJ 
Corollary[2] and Lemma|4] in this section. The rate distribution 
is validated by sweeping over a range of rate thresholds. The 
rate distribution obtained through simulation and that from 
Theorem [TJ and Corollary [2] for two values for the pair of bias 
and resource partitioning fraction (B,ij) is shown in Fig. 2a 
The respective densities used are Ai = 1 BS/km 2 , A2 = 5 
BS/km 2 , and A u = 50 users/km 2 with ai = 3.5, 02 = 4. 
The assumed transmit powers are Pi = 46 dBm and P 2 = 
26 dBm. Thermal noise power is assumed to be a 2 = 10 
dB. The rate distribution for the case with limited backhaul 
obtained through simulation and that from Lemma [4] is shown 
in Fig. [2b] The rate distribution is shown for two different 
backhaul bandwidths for a bias of B = 10 dB and without 
resource partitioning. Both the plots show that the analytical 
results, Theorem [T] and Lemma |4] give quite accurate (close 



to simulation) rate distribution. Furthermore, the mean load 
approximation based Corollary [2] is also not that far off from 
the exact curves in Fig. 2a This gives further confidence that 
the rate distribution obtained with mean load approximation 
in Corollary [2] can be used for further insights (as is done in 
the following sections). 

IV. Insights on optimal SINR and Rate Coverage 

As it was mentioned earlier the extent of resource parti- 
tioning and offloading needs to carefully chosen for optimal 
performance. Although a simplified setting is considered in the 
following results for analytical insights, it is shown that these 
insights extend to more general settings through numerical 
results. 

A. SINR Coverage: Trends and Discussion 

As stated before, the SINR coverage with and without 
resource partitioning exhibits different behavior in conjunction 
with offloading. The following lemma presents some key 
trends for SINR coverage in both settings. 

Corollary 3. Ignoring thermal noise (a 2 — > 0), assuming 
equal path loss coefficients and equal to four = 4), the 
SINR coverage without resource partitioning is 



1 



1 



V^tan-^Vr 7 ) + 1 + ^(^tan- 1 ^) + s/W) 



(36) 



where b = g 3 -, a = and p — p 3 -. The SINR coverage with 
resource partitioning for the corresponding setting is 

S(r) = ' 
+ TTtan- 1 ^) + 1 + ^(Vrtan-HVr) + 1 



N /rtan- 1 (Vr) + 1 + o^/p(Vr tan" 1 (^^76) + Vb) 
1 



/ftan-Hv^ + l + STp 



/?tan- 1 (^) + l+^ 



Proof: Using 
Q(t,4,x) = 1 + v^ 



du 



(37) 



= 1 + s/rt&n 1 (\/t/x), 



in Corollary |T| and substituting 



l , and b 



the expression in (37i is obtained. For the case with no 



resource partitioning, the expression derived in Lemma 5 of 



1 23 1 can be simplified using similar techniques to give ( 36 1 



Moreover, in this setting the following three claims can be 
made: 

Claim 1: Offloading with a bias (b > 1) leads to 
suboptimal SINR coverage S w in the case of no resource 
partitioning and t > 1 (0 dB). 

Claim 2: With resource partitioning, the bias b max- 
imizing the SINR coverage S can be greater than 
dB, the upper bound on which, however, decreases with 
increasing density of small cells. 
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Claim 3: With resource partitioning, the SINR coverage 
obtained by offloading all the users to small cells, i.e., 
b — » oo, is always less than that of no biasing, i.e., 6 = 1. 

Proof: See Appendix [D] ■ 
From the above corollary, it can be noted that S\b=i = 
S w \b=i, i.e., with no biasing/offloading SINR distribution with 
and without resource partitioning are equal as the orthogonal 
resource is not utilized by any user. Further, with resource 
partitioning the contribution to coverage from range expanded 



users (third term in (37i) increases with increasing bias, 
whereas the corresponding contribution from the macro cell 



users (first term in (37i) decreases with increasing bias. This 
is a bit counter-intuitive as one would expect with increasing 
bias only very good geometry users remain with macro cell. 
The reasoning behind this is the fast decrease in the fraction 
of such users A\ with increasing bias, which leads to an 
overall decrease in the coverage contribution from macro cell 
users. Similarly, the corresponding fast increase in the fraction 
of offloaded users Ab leads to an overall increase in the 
their contribution to coverage. A similar trend is observed for 
coverage without resource partitioning in ( [36] ). 

An intuitive explanation for the claims is as follows. Claim 1 
states that without resource partitioning, proactively offloading 
a user to small cell through an association bias is suboptimal, 
as the user would then always be associated to an AP offering 
lower SINR. On the other hand, with resource partitioning 
certain fraction of users can be offloaded to small cells and 
served on the resources which are protected from macro cell 
interference. In this case, increasing the small cell density, 
however, increases the interference on the orthogonal resources 
for the offloaded users, and hence, the optimal offloading bias 
is forced downward. Claim 2 justifies the described intuition 
when the bound is tight. Claim 3 suggests that preventing 
only offloaded users from macro tier interference is clearly 
suboptimal when almost all the macro cell users are offloaded 




8 10 12 
Association bias (dB) 



Fig. 3: Effect of small cell density on SINR coverage, with 
and without resource partitioning, as association bias is varied. 



to small cells. Of course, with large association bias, 6^oo, 
it would be better to shut the macro tier completely off, i.e., 
r) = 0. 

The above discussion is corroborated by the results in Fig. 
[3] which shows the effect of association bias on SINR coverage 
with varying density of small cells for a setting with ot\ = 3.5, 
a 2 = 4, and r = 0.5 (—3 dB). Without resource partitioning, 
it can be seen that any bias is suboptimal whereas for the case 
with resource partitioning optimal SINR coverage decreases 
with increasing density and the optimal bias also decreases. 
In case of no resource partitioning, increase in coverage with 
density is observed due to higher path loss coefficients of small 
cells (this was also observed in [18]). Thus, as evident from 
these results, the above claims hold in general settings too. 
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B. Rate Coverage: Trends and Discussion 

The following Corollary provides the rate coverage expres- 
sions for a simplified setting, which is used for drawing the 
following insights. 

Corollary 4. Ignoring thermal noise (a 2 — » 0), assuming 
equal path loss coefficients and equal to four (ctk = A), the 
rate coverage without resource partition and with the mean 
load approximation is 



TV 



+ 



1 

j~[ta,n~ 1 (^/u^) + 1 + a^/p{^/u{ta,n~ 1 (^/u 1 /b) + Vb) 
1 



(38) 



where Uk — t(pNk), b = g 2 , a = j 2 -, and p = p^. The rate 
coverage with resource partitioning under the corresponding 
assumptions is 



K = 



1 



mi tan 1 (^/uT) + 1 + a % /p(y^uT tan 1 (^m/b) + \fb) 
1 

M^tan-^^w^) + 1 + 5 ^p(tan _1 ( v %i) + 1) 
1 



Iub tan {*Jub) + 1 
1 



i 

a\/pb 



^/u^ tan 1 (0Tb) + 1 + ^ : 
where u\ = t(pNi~/i). 



(39) 



Proof: For the case with resource partitioning, the rate 
coverage expression follows from Corollary [2] using similar 
techniques as in the proof of Corollary[3] Without any resource 
partitioning Corollary 2 of J23~) is used. ■ 
From the above expressions it can be observed that 
7^1)7=1,6=1 = 72."' 1 6=1 > i- e -' the rate distribution for both sce- 
narios is same when no orthogonal resource is made available 
and there no offloaded users. For the case with no resource 
partitioning, the contribution to rate coverage from macro cell 



users (first term of (38 i) increases initially with increasing bias 
as the number of users sharing the radio resources at each 
macro BS decrease. But, beyond a certain association bias, 
due to the decreasing fraction of macro cell users, the overall 
contribution of the corresponding term towards rate coverage 
decreases. Similar trend is shown by the contribution from 
small cell users (second term in (|3~8|l). The initial increase with 
bias is due to the increasing fraction of small cell users and 
the subsequent decrease is due to increased number of users 
sharing the radio resources. This behavior of rate coverage 
could be seen as an intuitive reasoning behind the existence 
of an optimal bias. 

With resource partitioning, increasing rj increases the rate 
coverage of macro cell users and small cell users (first two 



terms of (39 1 respectively), whereas that of range expanded 
users decreases (last two terms of (|39|)), due to the decrease 
in available radio resources. With increasing bias, the rate 
coverage contribution from small cell users remains invariant 



(second term in ( 39 1) as the set Ug is independent of associ- 
ation bias. The contribution to rate coverage from the macro 
cell users (first term in ( f39] >) and that from offloaded users 
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Fig. 4: Effect of association bias, B, on rate coverage with 
A 2 = 5Ai. 



(sum of third and fourth term in ((39])) show similar variation 
with association bias as in the case of no resource partitioning. 
Therefore, there should exist an optimal bias for each T) in this 
setting too. 

The discussion in the above paragraphs is extended further 
in the following sections where the impact of various factors 
is studied on optimal offloading. For the following results, the 
parameters used are the same as in Section |III-E| with rate 
threshold p — 500 Kbps wherever applicable. 

1 ) Impact of resource partitioning: The effect of associa- 
tion bias and resource partitioning fraction on rate coverage, 
as obtained from Theorem [T| is shown in Fig. |4] The optimal 
pair for this setting is B = 15 dB and r\ = 0.53 (obtained by 
two-fold search), which gives a significant increase in the rate 
coverage compared to the case with no resource partitioning 
and offloading (B = dB, rj = 1). With decreasing resource 
partitioning fraction, the optimal association bias increases as 
more resources (macro interference free) become available for 
offloaded users. The variation of rate coverage with resource 
partitioning fraction for different association biases is shown 
in Fig. [5] The optimal resource partitioning fraction decreases 
with increase in association bias as more resources are needed 
to serve the increasing number of offloaded users. As shown, 
at lower association bias higher resource partitioning fraction 
is better as there are not enough range expanded users to take 
advantage of the resources obtained from muting the macro 
tier. 

A trend similar to rate coverage can also be seen in the 
5 th percentile rate p^ (where TZ(p^) = 0.95, i.e., fifth 
percentile of the population receives rate less than P95) in 
Fig. [6] The corresponding effect on median rate is shown 
in Fig. [7] The optimal pair of (B,7j) for these two metrics 
is same as that in rate coverage result. This shows that a 
single choice of the operating region provides a network-wide 
optimal performance across different metrics. 

2) Impact of infrastructure density: The impact of density 
of small cells on the fifth percentile rate is shown in Fig. [8] 
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Effect of resource partitioning fraction, rj, on rate Fig. 6: Effect of association bias and resource partitioning 
e with A2 = 5Ai. fraction (B,rj) on fifth percentile rate. 



It can be observed that at any particular association bias, as 
small cell density increases, ,095 also increases because of the 
decrease in load at each AP. With no resource partitioning, 
rj = 1, the optimal bias is seen to be invariant (at 5 dB) 
to the small cell density. Similar trend was also observed in 
p7) through exhaustive simulations. However, with resource 
partitioning, rj < 1, optimal association bias decreases with 
increasing small cell density. The optimal resource partitioning 
fractions (also shown for each density value) increase with 
increasing small cell density. These observations regarding the 
behavior of bias and resource partitioning fraction with small 
cell density can be explained by re-highlighting the learning 
from Sec. |IV-A| about optimal bias for SINR coverage. Without 
any resource partitioning, the optimal bias for SINR coverage 
is dB and independent of small cell density and similar 
independence is seen for rate coverage where the optimal bias 
is 5 dB. The insight of strictly suboptimal performance by 
a positive bias from SINR, though is clearly not valid for 
rate. With resource partitioning, increasing small cell density 
decreased the SINR coverage due to the increased interfer- 
ence in the orthogonal time/frequency resources allocated to 
range expanded users. Similar decrease of optimal association 
bias with increasing small cell density is seen for rate for 
same reasons. The increased interference in the orgthogonal 
resources also leads to the decrease in the optimal fraction of 
such resources, 1 — rj. 

It is worth pointing here that, although, the optimal asso- 
ciation bias decreases with increasing small cell density, the 
optimal traffic offload fraction A2 increases. With increasing 
density, at the same traffic offload fraction, the load per AP 
decreases, increasing the affinity of users for the corresponding 
tier. This trend is shown in Fig. [9] for two cases - (1) 
infinite backhaul bandwidth and (2) limited small cell backhaul 
bandwidth O2 = 5 Mbps. Decreasing backhaul bandwidth 
lowers both the optimal bias and consequently optimal offload- 
ing fraction at a fixed density. Interestingly, with increasing 
density, the optimal bias for case (2) increases, in contrast to 
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Fig. 7: Effect of association bias and resource partitioning 
fraction (B,rj) on median rate. 



the trend in case (1). This is because of the decreasing load 
per AP and the, thus, diminishing performance gap between 
the limited and infinite backhaul setting. Therefore, beyond a 
certain density the gap almost vanishes and the corresponding 
optimal biases reduce with increasing density. 

V. Conclusion 

We have provided an analytical framework for offloading 
and resource partitioning in two-tier co-channel networks. 
With the heterogeneity in base stations becoming prominent 
in emerging wireless networks these two techniques are set to 
play an important role in radio resource management. To the 
best of our knowledge, this is first work to derive rate distribu- 
tion in a heterogeneous cellular network, while incorporating 
resource partitioning and limited bandwidth backhauls. The 
availability of a functional form for rate as a function of system 
parameters opens a plethora of avenues to gain design insights. 
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Fig. 8: Variation in fifth percentile rate with association bias 
and resource partitioning fraction (B , 77) for different small 
cell densities. 
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Fig. 9: Effect of backhaul bandwidth and small cell density 
on the optimum association bias and optimum traffic offload 
fraction. 



Using the developed analysis, the importance of combining 
load balancing with resource partitioning was clearly estab- 
lished. It was further shown that the rate is a key metric for 
studying these techniques and insights based on just SINR are 
inconclusive. Investigating the impact of offloading on rate 
distribution in the uplink of HetNets could be an interesting 
future research direction. 
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p(Zi > (Pi) 1/ai z 1/& Afz 2 (z)dz 

>0 v ' 

(p 2 b 2 z 2 - q2 > Pizr Q1 n P2Z 2~ Q2 < p i z i ai 




1/ai 



Z 1/&1 <Z!< (Px 



fz 2 {z)dz. 
(40) 



Now 



P (Z k > z) = P ($ fe n 6(0, z) = 0) = cxp (-nX k z 2 ) , (41) 

where 6(0, z) is the Euclidean ball of radius z centered at 
origin. The probability distribution function (PDF) fz k (z) can 
then be written as 

fz k (z) = ^-{1 - ¥(Z k > z)} = 2n\ k z exp(-7rA fc z 2 ), Vz > 0. 

(42) 



Using (41 1 and (42i in (40 » gives Lemma [T] 



Appendix B 

Proof of Lemma ^ As a random user is more likely to 
lie in a larger association region then in a smaller association 
region, the distribution of the association area of the tagged 
AP of tier-fc, C k , is proportional to its area (under area 
approximation) and can be written as 

fci(c)ocefc h (c). (43) 

k 

Using the normalization property of the distribution function 



and (14 1 along with the distribution proposed in [28] for C(), 



the biased association area distribution is 

cfc k (c) 



E[C fc ] 
3.5 3 5 A fe 



A. 



3.5 



exp 



(44) 



r(3.5) A k \A k 

For a typical user u € Ui, the load at the tagged AP Ni 
comprises of the typical user and other users N a j (say) 
associated with the AP. The location of the other users in the 
association region of the tagged AP follows the reduced Palm 
distribution of $„ which is the same as the original distribution 



since <P U is a PPP |25 Sec. 4.4]. Using the properties of a PPP 
(30) and (g4j, the PMF of load is 

(X c )(™-i) 

exp(A„c)— — -/ c -(c)dc 

00 \ n L J- 1 



Pl (n)=W(Ni = n) 
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T(n + 3.5) (\ u Ai 



(n-1)! T(3.5) V A 



V(i) 



3.5 



^uAi 
A 



-(n+3.5) 



n>l. (45) 



Appendix A 

Proof of Lemma [TJ Using the definition of the three 
disjoint sets, the respective association probabilities are 



A x 



(PiZf" 1 > P 2 B 2 Z 2 - Q2 
P (z 2 > (P 2 B 2 

.n V 



f Zl (z)dz 



Appendix C 

Proof of Lemma [?} In this proof we first derive the 
distribution of the distance between the typical user u and 
the tagged AP when u €Ui. Let Y/ denote this distance, then 

( z j(i) > v,ug ui) 
¥{ueUi) 

(46) 



nY l >y)=¥(Z, m >y\ueUi) = 



12 



Using the proof of Lemma [T] in Appendix [A] 

;cxp -7ry^ Afc(PfcB fc j~ 



P(Yx > y) 
2ttAi 



F(Y B > y) 



2 expf-^^A fc (P fc ) 2/a ^ 2/<ife V 2 
Ab J *>y \ k=i J 



(47) 
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zexp \^TvJ2MPk) 2/ak z 2/&k jdz 

The above expressions lead to the corresponding PDFs 
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fvAv) 
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fc B (y) 



At 
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yexp -^^A fc (P fe B fe ) 2 /"^ 2 /^ 

V fe=i 

2/ex P (-^^A fe (P fe ) 2 /^y 2 / Sfc ) 



/jexp ^-Tr^A^PfeBfe) 2 /"^//"^ 



exp -^A fe P 2 / a V /afc 



fe=i 



(48) 

Now the moment generating function (MGF) of interference 
is derived. Expanding the interference, the MGF is given by 



M Ik (s) =E Ik [cxp(-s4) 
= E 



exp ( -sP fe H **~ ak 



(a) 



n m ^ 



(sP k x~ ak ) 



(49) 
(50) 

(51) 



Using change of variables with t = (sPfc) 2 / ak x 2 , the integral 
can be simplified as 



2x 



1 + (sP fc ) 



dx = (sP k ) 2/a " f 



At 



-2/c kz 2 1 + t a k/ 2 



(sP k ) 2/a "Z(l,a k , 



sP* 



giving the MGF of interference 



/fc ( S )=exp(-^A fc ( S P fc ) 2/Qfc Z( 1, 



where 



POO 

Z(a,b,c) = a 2 / b / - 



di/ 



sPfc 



(55) 



, (56) 



The SINR coverage of user u € is 

5,(r) = / P(SINR(y) > r|u € W,)/y, (y)dy. (57) 
Using the SINR expression in Q 

(SINR(y) > t\u e Ux) = P I ^ 2 yy ->t\ (58) 



H y >y^P 1 - 1 T\j2l k +a 2 



= E 

(a) 



cxp^-y^rPr 1 |E / fe + a2 |) 



(59) 



(60) 



exp 



SNRi(y) 



f[E Ik [exp^y^rPr 1 ^)] ( 61 ) 



fe=i 

2 



^^'-SMR^jllMU^Pr 1 ), (62) 

where SNRi(y) = Piy 2 1 and (a) follows from the indepen- 
dence of Ik- Similarly 



= exp ^-27rA fc jf {l - M Hx (sP fe :r- Qfc ) } zda^ P(SINR(y) > r|« G W s ) 



(52) 

-P(- ^ / l + ( s p~)-i^ (l ' J ■ < 53 > 



exp ( ) TT M /s (y^rPr 1 ) , (63) 



and 



where (a) follows from the independence of if^, (b) is obtained 
using the PGFL (25) of $ fc , and (c) follows by using the MGF 
of an exponential RV with unit mean. In the above expression, 
Zk is the lower bound on distance of the closest interferer in 
fc* tier which can be obtained by using {I} as 



P(SINR(y) > t\u e U B ) 



exp 



SNR 2 (y) 



/ 2 (y" 2 ^ 1 ) • (64) 



if 1 = 1 :z 2 = {P 2 B2) 1/a2 y ai/a2 ,z 1 = 
if; = S:z 1 = (P 1 ) 1 / Q2 y Ql / Q2 ,z 2 = y 
if I = B : Z2 = y 



(54) 



Using the PDF distribution (gSJl in (|57j along with (|62|)-(|64T) 
and ( |56| , the SINR coverage expressions given in Lemma[3] 
are obtained. The overall SINR coverage of a typical user is 
then obtained using the law of total probability to get S(t) = 
Ei<Si(r)A t . " ■ 
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Appendix D 

Proof of claims: Claim 1: The partial derivative of S u 
with respect to offloading bias b is 



-ay/p 



- — — + — 

b+T 2Vb 2Vb 



l+rb 2b 3 / 2 



{^tan- 1 ^) + 1 + ^(v^tan-^Vftf) + y/l/b)}' 



Since b.r > 0, hence -j^-^ + -K? > 0. Also, if 



t > 1 t > 1/6 for b > 1 
r 1 

> = 



T > 



1 



(65) 



1 + Tb 26 3 / 2 



26 3 / 2 - 6 
V b 5 lu < 0. (66) 



Thus, for t > 1 the SINR coverage decreases for all 6 > 1. 
Claim 2: Approximating tan -1 (a) s» a and substituting x for 
\/6, the partial derivative of coverage with respect to x is 

V B 5=VxJ 



+ ^tan- 1 (^) + l + ^} 
<Vp(£~1) , 1 



(67) 
(68) 



where v = -^/rtan 1 (\/t) + 1. The roots of the equation 
V x iS = are the zeros of the polynomial 



P(x) = x A a 2 p{v A - 1) + 2a: 3 a v /pw(l - r) 



c 2 {v 2 — 1 + a 2 pT(v 2 + 2)} — 4a;ay / puT — a 2 pr 2 — r. 

(69) 

Since u > 1, using the Descartes sign rule the polynomial 
P{x) has 1 positive root and upto 3 negative roots. The value 
of the positive root can be be upper bounded pT) by 
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1/4 



1/2 



ifr>l. (71) 



Further, the upper bound on the positive roots of P{—x) is 
given by 



L = max 



v 2 - l + a 2 pr(v 2 + 2) 



-,1/2 



a 2 p{v 2 — 1) 



2 2 

a pr + r 



1/4 



if r > 1 (72) 



20^/^(1 - r) 
a?p(v 2 — 1) 



d 2 - l + a 2 pr(u 2 + 2) 



-,1/2 



2 2 i 

a pr + r 



a?p(v 2 — 1) 
1 1/4 



T?p(v 2 — 1) 



ifr<l. (73) 



Note that — L is the lower bound on the negative roots of P(x), 
since they are same as the positive roots of P(—x). Clearly, 
both U and L are inversely proportional to the density of small 
cells a. Since 



-L<Vb<U =^ b < max{[/ 2 ;J L 2 } , 



(74) 



therefore the upper bound on optimal bias is inversely propor- 
tional to the density of small cells a. 

Claim 3: The SINR coverage at very large offloading bias is 



51, 



= lim S 

b—>oo 

i i 



V V + 



(75) 
(76) 



where v = yVtan 1 (y/r) + 1. With the knowledge that 

5| fc= i = S w \ b= x we get 



S\b=l — S\b=oo — 



1 



1 



1 



Thus, 51 



<<S! 



v + av v w+i v(v + i)(« + av) 
> since v > 1. 

b=i. ■ 
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